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Abstract 

We study transport properties of a strongly correlated monoatomic chain coupled to metallic leads. Our system 
is described by tight binding Hubbard-like model in the limit of strong on-site electron-electron interactions in the 
wire. The equation of motion technique in the slave boson representation has been applied to obtain analytical 
and numerical results. Calculated linear conductance of the system shows oscillatory behavior as a function of the 
wire length. We have also found similar oscillations of the electron charge in the system. Moreover our results show 
spontaneous spin polarization in the wire. Finally, we compare our results with those for non-interacting chain 
and discuss their modifications due to the Coulomb interactions in the system. 
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1. Introduction 

Recently one-dimensional (ID) quantum wires 
(QW) have attracted much attention due to their 
potential applications in nanoelectronics [1] and 
quantum computing [2]. The knowledge of the 
transport properties of such structures is crucial 
for the design and fabrication of the nanodevices. 
On the other hand the quantum wires, although 
conceptually simple, are very interesting from sci- 
entific point of view as they display extremely rich 
phenomena, very often different from those in two 
and three dimensions [3,4]. The understanding 
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of the properties of such ID objects is a major 
challenge in the field of nanophysics. 

The conductance of the quantum wires has been 
studied both experimentally and theoretically by 
number of authors (see [5] for a review) . The exper- 
imental studies require advanced techniques of fab- 
rication of such structures. Those include: growing 
of QW on metallic surfaces [6,7], scanning tunnel- 
ing microscope techniques [8] or mechanically con- 
trolled break junctions [5,9,10]. Those fabrication 
techniques allowed for revealing of many phenom- 
ena like charge quantization in units of Go = 2e^//i 
[11], deviations from that (0.7(2e^//i) anomaly) 
[12], spin-charge separation (Luttinger liquid) [13], 
oscillations of the conductance as a function of the 
length of the chain [9,10] or spontaneous spin po- 
larization in QW [12,14]. 

The purpose of the present paper is two fold. 



Preprint submitted to Elsevier Seience 



2 February 2008 



The first one is to investigate the oscillations of the 

conductance as a function of the wire length in the 
case of strong Coulomb interactions. The oscilla- 
tory behavior of the conductance manifests itself 
as a maximum of the conductance when a number 
of the atoms in a wire is odd and minimum when 
the number is even. This effect is known as the 
even-odd conductance oscillations. Most common 
examples are the oscillations with a period of two 
[9,10,15,16,17,18] and four atoms [19]. However, 
the conductance can oscillate with different (from 
two and four) periods, depending on the average 
occiipation of the wire. Moreover, recently the an- 
alytical formulas for M-atom (M > 2) oscillations 
have been found [20]. However, those analytical 
formulas remain valid for noninteracting wire only. 
In the presence of strong Coulomb interactions the 
even-odd oscillations with a period of two atoms 
have been also found [21,22,23,24,25,26]. The M > 
2 conductance oscillations have been only reported 
for the nearest neighbor Coulomb interactions [27] . 
Therefore we shall study the oscillations of the con- 
ductance in the case of strong on-site Coulomb in- 
teractions and see how their period will be modi- 
fied. 

The second purpose is to see if the wire will 
exhibit any spontaneous spin polarization in the 
presence of strong correlations. Such spontaneous 
polarization has been observed experimentally 
[12,14]. But it is well known that the ferromag- 
nctism in strictly ID objects is forbidden due to 
the Lieb-Mattis theorem [28]. 

The paper is organized as follows. In Sec. 2 we 
present theoretical description of the model wire, 
in Sec. 3 we show the results of the calculations 
of the conductance, charge and spin polarization. 
Finally, in Sec. 4 we provide some conclusions. 



2. Theoretical description 

Our system consists of the quantum wire mod- 
eled as a chain of N atoms coupled to the left L 
and right R lead described by the following Hamil- 
tonian in the limit of strong on-site Coulomb in- 
teraction (Ui — » do) in the slave boson represen- 
tation where the real wire electron die is replaced 



by the product of the boson bi and the fermion fi^^ 
operators {dia = bf fi^) [29,30,31]: 

Ak(7 id 

+ E ^L{R)i^(^t{R)kaK(N)fHN)a + h.C, (1) 
keL{R),cT 

where caicct stands for the electron with the single 
particle energy e^kj the wave vector k and the spin 
a in the lead A = L,R. Ei denotes the wire en- 
ergy level at site tij is the hopping integral of the 
electrons between neighboring wire sites i and j, 
and Vi;,(i{)k is the hybridization matrix element be- 
tween electrons at site 1(A'') and those in the lead 
L{R). 

In the linear response and at the zero tempera- 
ture the conductance G is proportional to the total 
transmittance T, i. e. G = ^T^- our case the 
transmittance is given by [20,32] : 

Tm{E)=Y,^lTr\G\^,{E)\\ (2) 

cr 

where G\^^ is the retarded Green fimction 
(GF) connecting the ends of the wire (sites 1 
and A'') and ^l(r) is the elastic rate ^l{r) = 
^T^Sk \^L{R)Vi?S(E - enB)v). In calculations we 
have assumed constant bare density of states in 
the leads. 

Using the equation of motion technique for the 
retarded GF with Hamiltonian (1) one can write 
the general matrix equation for G^^^ in the form: 

I^g; - N^. (3) 

Due to the strong on site Coulomb interactions in 
the wire Ui, which is assumed to be infinity in our 
case, the problem cannot be solved exactly and 
one has to make approximations of the higher or- 
der GFs emerging in the equation of motion for 
the retarded GFG[^-, = ((&+/,,|/+ We have 
used Hubbard I like approximation [31] accord- 
ing to which the GF ((/+ 6,))^ is 

approximated by {f^_eft-a){{bk fka\f^„bj))E and 
the other higher order GFs are neglected. This ap- 
proximation is reasonably good for not very large 
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values of the hopping t and neglects higher order 
processes, like for example the Kondo effect. 

Within the present approximation scheme Aa- in 
the Eq.(3) is N x N tridiagonal symmetric matrix 
with the elements: 

r 

Act = {E - Si)6ij + i — {6i,l6j,l + 6i,NSj,N) - 

t[{l - ni-^)6i^i+i + (1 - ni+i-a)5i+i^i] (4) 
and Na is the diagonal matrix of the form: 

Ncr = {I - ni^a)Stj , (5) 

with ma = {ftM = -I j dEImGlaiE) being 
the average occupation of the electrons with spin 
a at site i. 



3. Results 

In numerical calculations we have assumed all 
the wire site energies to be equal (sj = eq) and 
similarly hopping integrals tij = t. All energies are 
measured with respect to the leads Fermi energy 
Ep = in units of F = Fi, = F/j = 1. Moreover, 
the occupation is calculated self-consistently 
on each wire site. 

To find the condition for M-atom conductance 
oscillations one has to solve the relation: T/v = 
Tn+m, where Tjv {Tjsi+m) is the transmittance of 
the wire consisted oi N {N + M) atoms, given by 
Eq. (2). In general, for t/j ^ oo it is not possible to 
get analytical expression for the oscillations con- 
dition without further assumptions. Note that for 
Ui = Q the problem can be solved exactly and such 
condition can be found [20]. In this case it reads 
•^o-^'dr) = where Z = 1, 2, ... M - 1. 

In the case of [/, ^ cxo both matrices A^ and 
Na depend on the occupation on each wire site 
rii^a- thus the problem has to be solved numerically. 
However, if one assumes that the occupation is the 
same on each wire site and does not depend on spin, 
i. e. riia = ria = n-a, similar M-atom oscillations 
condition as for Ui = can be found. In this case 
it reads: 



with I = 1, 2, ... M — 1. Unfortunately, one has 

to know the average occupation ria ■ The only case 
of M = 2 can be solved analytically. As one can 
see from Eq. (6) the period of two can be obtained 
for Ep — eo = 0, i.e. when the wire single particle 
energy levels all coincide with the Fermi energy. 

In Fig. 1 the total linear conductance G — 
X^cr is plotted as a function of the wire length 
N and the energy level eo for t — 4. As one can 
see all the figures show similar patterns with the 
regions of large and small conductances. More- 
over, for £0 = behavior of the conductance is 
the same in all three cases which leads to the 
conclusion that correlations arc not important in 
this case. It always shows even-odd {M = 2) os- 
cillations. Away from eq = correlations strongly 
modify the conductance, shifting the maxima of 
G (for fixed N) towards lower (higher) energies 
for Eo > (eo < 0). This is due to the modifi- 
cation of the wire hopping, which depends now 
on the occupation rua (see Eq. (4)). This effect 
also leads to the strong asymmetry for positive 
and negative energies. For negative (positive) wire 
energies the occupation is large (small), therefore 
effective hopping i = t{l — rii-a) is small (large), 
thus the conductance decreases (increases). Note 
that there is no such asymmetry in the case of 
Ui = 0. Moreover, different periods of the conduc- 
tance oscillations can be observed, depending on 
the position of the wire energy level. 

Another important finding is that the wire shows 
spontaneous spin polarization. It is well known 
that in strictly ID wire the spin polarization is 
prohibited due to the Lieb-Mattis theorem [28]. 
However this theorem is valid for infinite wire only. 
In experimental situation the wire is always con- 
nected to the electrodes and this is why the spin 
polarization is observed experimentally [12,14]. In- 
terestingly, it was predicted recently that even infi- 
nite wire, but of the zig-zag shape, can also exhibit 
the spin polarization [33] . 

The total conductance in ferromagnetic case is 
shown in the bottom panel of Fig. 1. Again, unlike 
for J7, = 0, the conductance pattern shows strong 
asymmetry for positive and negative energies. Note 
that for negative energies there are no differences 
between G/erro (bottom panel) and Gpara (mid- 
dle panel). For such energies iterations always con- 
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to spend more time on the same site than move to 
another one. Moreover in the case of the lack of the 
inter-site interactions any collective phenomenon 
is not possible. Situation is different for positive en- 
ergies. In this case the electrons are more mobile, 
as the effective hopping is larger due to the small 
values of the wire occupations. Thus they interact 
with each other via hopping and it is possible and 
energetically more favorable to get the ferromag- 
netic state. The hopping, which is the correction to 
the position of the wire energy levels, leads to the 
effective splitting of this levels. Thus the conduc- 
tance is spin polarized in this case. It can be read off 
from Fig. 2, where the difference between spin up 
and spin down conductance is displayed. It turns 
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Fig. 2. The difference between spin up and spin down 
conductance as a function of the wire length A'^ and the 
wire energy level eg- 



Fig. 1. The total linear conductance G = ^^^Gcr as a 
function of the wire length N and the wire energy level eo 
for Ui = (top panel), f/i — ► oo in paramagnetic configu- 
ration riio- = rii-cr (middle panel) and ferromagnetic one 
(bottom panel). 



verge to paramagnetic solution. This can be again 
explained by effect of the hopping modification. 
For negative energies the wire occupation is large, 
thus the effective hopping t = t(l — rii^^) is small. 
One can imagine in this case that the electrons are 
more localized and it is more convenient for them 



out that the strongest differences between G| and 
Gi can be found for intermediate values of eq- For 
energies close to Ep the modifications are weak due 
to the small values of the hopping while for higher 
energies the wire occupation is very small thus the 
difference between spin dependent effective hop- 
ping can be neglected. Finally it is worthwhile to 
note that for = 1 (single atom) there is no fer- 
romagnetic solution as the ferromagnetism is gov- 
erned by inter-site hopping in our case. Whether 
this spontaneous spin polarization is a true effect 
or a drawback of the approximation used remains 
an open question, as it is known that the mean 
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field like theories overestimate the role of the mag- 
netism. The problem will be further studied. 

Corresponding spin polarization — (n^ = 
J2i ''^i'j/N) is shown in Fig. 3. As one can see the 
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Fig. 3. The spin polarization as a function of the wire 
length A'' and the wire energy level eo- 

spin polarization pattern is similar to that of the 
conductance differences (see Fig. 2). 

In Fig. 4 we show the conductance (left panels) 
and the occupation (right panels) as a function of 
the wire length for a number of the energy levels eo 
in the paramagnetic configuration. The values of 
£o have been chosen in such a way that they lead 
to the maxima of the conductance for = 1, 2 
and 3 atom wire. For example, the maxima of the 
conductance of the N = 2 atom wire correspond to 
Eq ~ —2.298 and 3.135 (see the middle panel of Fig. 
1). These maxima in [/^ = case give the periods 
of the conductance oscillations. To find A/-atom 
period it is enough to determine the maxima of the 
conductance for iV = M — 1 atom wire [20] . As one 
can read off from Fig. 4, depending on Eq, one gets 
different periods of the conductance oscillations. 
Moreover, except for the special case of M = 2, 
the amplitude of the oscillations decreases with the 
wire length. This is a kind of damped oscillations. 
No such effect has been observed for Ui = [20] . 

At this point we would like to comment on the 
other results known in the literature. The even-odd 
(M = 2) oscillations problem was extensively stud- 
ied within the second order perturbation theory 




Fig. 4. The total conductance (left panels) and the occu- 
pation (right panels) vs. wire length. The left panels show 
the conductance with different oscillations periods (2 to 4) 
from top to bottom. The positions of the wire energy levels 
are indicated in the figure. 

in Ui (SOFT) [21,22] and the numerical renormal- 
ization group (NRG) approach [23,24]. The results 
show similar behavior of the conductance for odd 
number of atoms in a wire: it always reaches the 
unitary limit (^), independently of Ui. Such be- 
havior is a consequence of the Kondo effect. How- 
ever, in our case the situation is slightly different, 
as we get M = 2 oscillations in the mixed valence 
regime only (eq = 0) , where the Kondo effect is ex- 
cluded. Thus our even-odd oscillations are caused 
by the resonances associated with the energy level 
structure of the chain rather than the Kondo ef- 
fect. On the other hand, the conductance for even 
number of atoms in a wire is strongly suppressed, 
in agreement with SOFT [21,22] and NRG ap- 
proaches [23,24]. However, NRG calculations show 
that the conductance exponentially depends on Ui, 
and in the limit of Ui —> oo vanishes, contrary to 
our results, as we get non-zero values of G (see Fig. 
4), depending on the hopping t. The conductance 
vanishes in the limit of very large or very small 
values of t. Interestingly, when t = T/2, the con- 
ductance reaches the unitary limit and shows no 
oscillations, i.e. is equal to ^ for even and odd N. 

Figure 4 (right panels) shows the wire length de- 
pendent occupation which also oscillates with the 
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same period as the conductance does, except for 
the special case of Eq = {M = 2), where the occu- 
pation remains constant. Moreover, the Coulomb 
interactions Ui lead to the reduction of the occu- 
pation oscillation amplitude. Similar effect, albeit 
for small Ui, has been found within self-consistent 
Hartrec-Fock approximation [34] . 

In ferromagnetic case the situation is more com- 
plex. For positive values of Eq due to the splitting 
of the conductance maxima (see Fig. 1) no regu- 
lar oscillations have been observed. On the other 
hand; for co < one gets such oscillations but 
in this case the solutions remain always paramag- 
netic. Figure 5 shows the comparison of the con- 
ductance and the occupations in paramagnetic and 
ferromagnetic configurations for the wire consisted 
of five atoms. For negative energies, as discussed 
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Fig. 5. The total conductance in paramagnetic (left top 
panel) and in ferromagnetic (right top panel) configuration 
and corresponding occupations (bottom panels) vs. eo of 
five atom wire. 

before, there are paramagnetic solutions only. For 
N > I and eq > the ferromagnetic solutions 
emerge in certain energy regimes. In this case n-f ^ 
ni and resulting conductance peaks are split. In- 
terestingly, the spin polarizations occurs only in 
the regimes where the occupation has a large slope 
or the total conductance has a maximum. More- 
over, the splitting of the conductance leads to the 
fact that Gferro shows more maxima than Gpara- 
A number of maxima of Gferro is related to the 
wire length (number of atoms - N) and for odd N 
it gives (3A^ — l)/2 maxima, while for even N one 



observes 3A^/2 maxima. In paramagnetic case it is 
always equal to a number of atoms N in the wire 
(compare top panels of Fig. 5). 

Finally we would like to comment on the valid- 
ity of our approach. The present calculations com- 
pletely neglect the Kondo effect which is important 
at low temperatures and Ed < 0. We expect some 
modifications in this regime, as this effect leads to 
the corrections of the conductance of the order of 

. Thus our results apply for temperatures higher 
than the Kondo temperatme. On the other hand, 
we do not expect any qualitative modifications in 
the mixed valence and the empty regimes {sd > 0). 



4. Conclusions 

In summary we have studied the conductance os- 
cillations of the strongly interacting wire as a func- 
tion of the wire length. We have found that strong 
Coulomb interactions significantly modify the pe- 
riods of the oscillations showing strong asymmetry 
for negative and positive wire energy levels. They 
also lead to the suppression of the conductance 
with increasing wire length. There are no such ef- 
fects for noninteracting wire. Moreover, strong in- 
teractions lead to the spontaneous spin polariza- 
tion for positive wire energies, observed in experi- 
ments. 

Acknowledgements 

This work has been supported by the grant no.l 
P03B 004 28 of the Pohsh Committee of Scientific 
Research. T. K. thanks the Foundation for Polish 
Science for a Fellowship for Young Scientists. 



References 

[1] D. R. Bowler, J. Phys. Cond. Matter 16 (2004) R721. 

[2] A. Bcrtoni, S. Reggiani, Semicond. Sci. Technol. 19 

(2004) S113. 

[3] J. M. Luttinger, J. Math. Phys. 4 (1963) 1154. 

[4] F. D. M. Haldane, J. Phys. C: Solid State Phys. 14 

(1981) 2585. 

[5] N. Agrait et al, Physics Rep. 377 (2003) 81. 
[6] M. Jalochowski et al, Surf. Sci. 375 (1997) 203. 



6 



[7] M. Krawiec et ai, phys. stat. sol. (b) 242 (2005) 332. 

[8] A. Yazdani et ai, Science 272 (1996) 1921. 

[9] C. J. MuUer et ai, Phys. Rev. Lett. 69 (1992) 140. 

[10] R. H. M. Smit et ai, Phys. Rev. Lett. 87 (2001) 
266102. 

[11] B. J. Wees et ai, Phys. Rev. Lett. 60 (1988) 848. 

[12] K. J. Thomas et ai, Phys. Rev. Lett. 77 (1996) 135. 

[13] O. M. Auslaender et al, Science 308 (2005) 88. 

[14] B. E. Kane et ai, Appl. Phys. Lett. 72 (1998) 3506. 

[15] H. -S. Sim at al., Phys. Rev. Lett. 87 (2001) 096803. 

[16] N. D. Lang, Ph. Avouris, Phys. Rev. Lett. 81 (1998) 
3515; Phys. Rev. Lett. 84 (2000) 358. 

[17] E. G. Emberly, G. Kirczenov, Phys. Rev. B60 (1999) 
6028. 

[18] R. Gutierrez et ai, Acta Phys. PoL 32 (2001) 443. 

[19] K. S. Thygesen, K. W. Jacobsen, Phys. Rev. Lett. 91 
(2003) 146801. 

[20] T. Kwapinski, J. Phys.: Condens. Matter 17, 5849 
(2005). 

[21] A. Oguri, Phys. Rev. B59, 12240 (1999); Physica 
B284-288, 1932 (2000); Phys. Rev. B63, 115305 (2001); 
J. Phys. Soc. Jpn. 70, 2666 (2001). 

[22] Y. Tanaka, A. Oguri, J. Phys. Soc. Jpn. 73, 163 (2004). 

[23] A. Oguri, A. C. Hewson, J. Phys. Soc. Jpn. 74, 988 
(2005). 

[24] A. Oguri et al., J. Phys. Soc. Jpn. 74, 1554 (2005). 

[25] R. A. Mohna et ai, Phys. Rev. B67 (2003) 235306. 

[26] V. Meden, U. SchoUwock, Phys. Rev. B67 (2003) 
193303. 

[27] R. A. Mohna et al, Europhys. Lett. 67 (2004) 96. 

[28] E. Lieb, D. Mattis, Phys. Rev. 125 (1962) 164. 

[29] P. Coleman, Phys. Rev. B29 (1984) 3035. 

[30] J. C. Le Guillou, E. Ragoucy, Phys. Rev. B52 (1995) 
2403. 

[31] M. Krawiec, K. I. Wysokihski, Phys. Rev. B59 (1999) 
9500. 

[32] S. Datta, Electronic transport in mesoscopic systems, 
Cambridge Univ. Press, Cambridge (1995). 

[33] A. D. Klironomos et al, cond-mat/0507387' 

[34] T. Kostyrko, B. R. Bulka, Phys. Rev. B67 (2003) 
205331. 



7 



